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Abstract: We reconsider analysis of data on the cosmic microwave back- 
' ground on the largest angular scales. Temperature multipoles of any order 

. factor naturally into a direct product of axial quantities and cosets. Striking 

coincidences exist among the axes associated with the dipole, quadrupole, and 
Q^. octupole CMB moments. These axes also coincide well with two other axes 

independently determined from polarizations at radio and optical frequencies 
propagating on cosmological scales. The five coincident axes indicate physical 
correlation and anisotropic properties of the cosmic medium not predicted by 
. the conventional Big Bang scenario. We consider various mechanisms, includ- 

ing foreground corrections, as candidates for the observed correlations. We also 
consider whether the propagation anomalies may be a signal of "dark energy" in 
' the form of a condensed background field. Perhaps light propagation will prove 

^ . to be an effective way to look for the effects of dark energy. 

1 Introduction 

The cosmic microwave background is full of information. It is conventionally 
analyzed by expanding temperature correlations in spherical harmonics. Cur- 
rent literature is dominated by attention to I ^ 200. The region of I ~ 200 is 
predicted by models of causal correlation at decoupling. Multipoles of small^ 
"spin" I contain information on scales too large to possibly be physically cor- 
related in the standard Big Bang model. Yet it is well known that small I 

^Spin of the quantum type cannot be confused, so we can use the simple term for brevity 
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quantities do not obey expectations, There has been a burst of workpP suggest- 
ing reinterpretation, usually with introduction of new parameters, in order to 
verify the Big Bang. 

We maintain that the small-/ components play a scientifically pivotal role. 
They have the potential to falsify current hypotheses. Falsification is more 
powerful than verification, so it is interesting to re-investigate the low moments 
with an open mind. We first set aside correlations, and use the temperature 
map AT{0, (f>)/T directly, because it has far more information. 

The Z = 1, m = 0, ±1 or dipole mode is commonly attributed to our motion 
relative to the CMB rest frame. This is a fine hypothesis, but we wish to think 
afresh about it. We find it sobering to realize that the hypothesis has never been 
tested. Everyone expects some component of dipole from a boost of our system, 
and yet the sum of any number of dipoles is a dipole in a combined direction. 
The annual motion of the Earth about the Sun is an observable relative to a fixed 
dipole, and has been measured. However the balance of our proper motion is not 
an independent CMB observable. Not unrelated, the confirmation of the "great 
attractor" of local gravitational attraction correlated with the dipole seems to 
have gotten complicated: it remains an unfinished project. The relative size 
of the dipole above an isotropic background might be argued to confirm the 
orthodox interpretation. Yet confirmation of three I = 1 variables with three 
free parameters is not a test. 

If the dipole were the only observable, perhaps discussion would be moot. 
Yet for 30 years the attempt to resolve the relation with higher-spin moments 
was frustrated by their unexpectedly small magnitude. The relationships be- 
tween the \ow-l moments are not the predicted ones to this day. It turns out 
that the \ow-l multipoles contain additional group-theoretic information on di- 
rectionality of the radiation background. The dipole is simple because it is 
summarized by a magnitude and a direction in space towards Virgo. Presum- 
ably this was simply a random accident. It is quite surprising that the various 
directions associated with higher multipoles show mutual correlations in direc- 
tion which are not predicted by the basic Big Bang. This causes us to rethink 
the meaning of the low / multipoles, preferring to work from the bottom-up as 
data guides us, and without deciding here what the origin of the effects may be. 

Here is the plan of the paper. In the next Section we develop a decom- 
position of the multipoles based on symmetry analysis. There is a notion of 
invariantly preferred frames for any spin-/. The axial quantities obtained co- 
incide with those calculated independently by a numerical method. We also 
review other possibilities, based on projecting arbitrary products of angular 
momenta representations into smaller dimension. Section 3 reviews previously 
published work and compares it. We separate statistical analysis into Subsec- 
tion 3.4. This analysis is deliberately incomplete, and conservative, because the 
statistical significance of data depends upon the hypothesis. Section 4 discusses 
possible mechanisms for the new effects found. The most ready explanation 
would perhaps be local effects on the CMB, but the importance of CMB data 
for cosmological deductions makes local perturbations hard to accept. Moreover 
the existence of anisotropic correlations at optical frequencies, and independent 
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Figure 1: AithofF-Hammer projections of the front hemisphere — 7r/2 < I < -k /2 
of the CMB temperature map in the spin I = 2 quadrupole (left) and / — 3 
octupole (right) modes. The back side is trivial from the symmetry of the 
modes. 

radio frequency data, is not explained by any known model. We find a loophole 
in current reasoning related to dark energy. We suggest that propagation and 
polarization anomalies may be an excellent probe of dark energy as compared to 
the more familiar Hubble flow probe of the energy-momentum tensor. Perhaps 
light will prove to be an effective tool to look for the effects of dark energy. 
Concluding remarks are given in Section 5. 

2 Factoring Modes 

Let > be the projection of the quantity ^{0, (j)) into angular momentum I. 
A moment is defined by ipim —< l,rn, \ip{l) >, where \ l, m > are eigenstates 
of the angular momentum operators J^, Jz in the spin-Z representation. The 
temperature measurements are real numbers, which produces a constraint on 
the matrix elements: 

i^lm = i-ir^Lm- (1) 

What are the natural ways to derive an axial quantity from a higher multi- 
pole? A simple assignment of a vector quantity is the expectation 

<m\j\m>- (2) 

This proposal has two flaws: (1) it is quadratic, not linear in the data, and of 
doubtful physical interpretation. The quantity in Eq. [21is a map from represen- 
tations V'*(0 ~' spinl with dimensions of temperature-squared. In the 
product of two representations many new representations of high spin are cre- 
ated and vectors of temperature-squared have no immediate meaning. Perhaps 
for quadratic correlations the quadratic products of temperatures giving spin-1 
might mean something. But we are discussing the temperature. (2) The quan- 
tity is zero. The expectation value in Eq. |2] vanishes for real representations. 
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2.1 A Linear Map 



Here we teach ourselves how to factor ?-polc expansions. The problem of fac- 
toring is geometrical and would be much the same for any angular quantity. 

If the state |^(Z) > were an eigenstate of with eigenvalue m' in the frame 
where is diagonal, we could obtain a relation 

Here |m >= \l, m > is a concise notation. Correspondingly, if we rotated the 
operator J, 

J ^ J{a, P, 7) = Uiia, /?, j)JUlia, /?, 7), 

then we could determine the same relation for an eigenstate in an arbitrary 

frame related by Euler angles a, /3, 7. It follows that there is a linear map or 
"wave function" t/)^ which contains the information on frame-orientation: 

The ranges of indices is A; = 1 ... 3 and m — —l . . .1. We call this a wave function 
in analogy with the Bethe-Salpeter wave function, which is the matrix element 
of the operator creating a quantum state (analogous to J± = Jx ± iJy) and a 
fiducial vacuum state (analogous to |m = >.) We are not doing quantum 
mechanics and the physical observables exist at the level of this wave function, 
not the wave function-squared. Most importantly, <j)'^ is linear in the data and 
transforms like a vector in the index k and a rank I tensor in the index m. 



2.1.1 Completeness 

The map to 0^ is simply a coordinate transformation. Start with =< 
m\il) >, of course for fixed I. Make 

'^m ~ ^mm,''^rn''i 

^'^ - ^ <m\j'' \m' > . (4) 



■ mm 



Repeated indices are summed. Since (/)^ has two indices, and ipm' has one index, 

the transformation is not obviously invcrtible. But the assignment of indices 
simply records the division of the space. The inverse of F^^, exists and is 
denoted r™™ . Construct the inverse directly: 

r^'""r^„„ = 777^ < m| J'= \m" >< m"\ \m' >, 
= ^<m|J^|m'>, 

— ^mm' ■ (5) 
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Technically r™™ is a pscudoinverse (Moorc-Pcnrosc inverse) for the two indices 
m, k lying on the subspace of fixed I. The map to (j) is unitary for each I. 

2.1.2 Invariantly Preferred Frames 

Although (j)^ is neither symmetric nor Hermitian in general, it has an invariantly 
preferred coordinate frame - in fact, two frames. The frames are found by the 
singular value decomposition (svd) 

<t>t = f2 • (6) 

a 

The matrix elements and are "unitary" (i.e. orthogonal when real): 

^<«'r =5„„.; (7) 

m 

Y,eUet'r=Sc.a'. (8) 

k 

It follows that in a coordinate system rotated by e on index k and rotated^ 
independently by on index m the expansion of (p!^ (A^, A^, . . .A"*) and 
diagonal. This is what is meant by "an invariantly preferred coordinate system" . 
In the preferred coordinates, the frame Ska, meaning that (1, 0, 

. . . ), and so on. The rest of 0^ becomes a definition of the particular vector 
chosen to be a fiducial basis: the case (j)Q is a perfectly adequate "vacuum." So 
the meaning of the svd is that we can relate any state by a few invariants and 
a rotation to an invariantly determined reference frames. 

Consistency requires that the number of terms A does not exceed the di- 
mension of the smaller space. For I = 1 there is a degeneracy and the number of 
terms A = 2. In all cases ^ < 3. The matrices are not a complete "frame" 
but have the nature of a coset. The eigenvalues A" are invariant, real, and by 
conventions of phase, positive. The svd is the direct generalization to two inde- 
pendent spaces of the usual fact that a Hermitian matrix can be diagonalized 
with the conjugate coordinate transformations acting on the left and the right. 
For those with backgroimds only in quark physics, the svd is the procedure to 
reduce the CKM matrix to its invariant terms and unitarily-related frames. 

Degeneracy If two or more eigenvalues A" are equal then tjj is invariant 
under the corresponding orthogonal symmetry. While this is rare, a degeneracy 
happens to occur in the data at hand, as seen momentarily. 

2.1.3 Axial Quantities 

Starting from our standard coordinate frame {\i >} with names i, with < j\i >= 
^ii) < j|3 >= (0, 0, 1), and so on, the three vectors e\, e^, e\ represent three 

^ "Rotation" is used figuratively because u lies on U{1) in the general case. We did not 
investigate conditions for the rotation on index m to literally lie in an SU{2) rotation subgroup. 
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preferred orthogonal axes in space. Sometimes we will use the notation for 
this reason. We assess the importance of these axes by the weights A". They 
parametrize the degree of asymmetry. Ordering the eigenvalues from largest to 
smallest, the first frame element designated by a = 1 describes a preferred axis. 
We reiterate that eigenvectors describe an axis, in the meaning of an oriented 
line in space, not signed-vectors in the sense of having a direction. In detail 
this comes because the signs of the e" are conventions fixed by the signs of the 
weights A". Numerical code can conceal this fact. 

2.1.4 The Quadrupole 

The quadrupole is very instructive. A quadrupole set of 5 numbers (9-2, 9-1, ■■■Q2) 
is just the relisting in the angular momentum basis of the 5 elements of a trace- 
less symmetric 3x3 Cartesian tensor Qij : 

= vr;^m{i). (9) 

Here V™ is a unitary set of coefficients to change indices m to indices ij. 
Such a tensor has an obvious preferred frame in which it is diagonal: 

3 

a 

How are the obvious frame elements > related to the frame elements |e" > 
defined by svd7 

They are the same frame, as evident by symmetry of diagonalizing on the left. 
The demonstration is awkward and we simply verified the fact with numerical 
examples. 

2.2 Frames Defined by Expectation Value 

De Olviera et al (dOZTH)[2| have obtained a special axis for the CMB quadrupole 
and octupole by the following procedure. Define n = n{9^ (/)) as a unit vector 
depending on angular parameters 6, (p. Find the maximum of the quadratic 
expectation of (J • n)^, 

<ij{l)\{J-nf\Tp{l)>~^max. (10) 

The sign of n is, of course, not determined by Eq. ^| dOZTH report finding 
maximal unit vectors by a numerical search with the 3,145,728 pixel map. Using 
galactic coordinates and the standard conventions of axes (see Appendix) , they 
find for the quadrupole and octupole: 

h{l = 2) = ±(-0.1145, -0.5265, 0.8424); 

n{l = 3) = ±(-0.2578, -0.4207, 0.8698). (11) 

We think this is an important result and we will return to it. 
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Besides one axis of absolute maximum weight, we observe there are always 
two other h per multipole giving local maxima. We see this by writing 

K^fe' =< J'^J'^' |V(0 > ■ (12) 

Since the commutator of the J's gives no contribution, K^^ is symmetric in 
its indices. By the Rayleigh-Ritz variational theorem fi will be one of three 
eigenvectors of K when the expectation is extremized. (There would be N 
such eigenvectors for N group generators.) That is, the variational problem 
determines a frame. 

Earlier we mentioned that tensor quantities quadratic in the temperature 
do not have the interpretation of temperature. What is the relation of the 
quadratically-produced K^^ axis and the one from svdl It is a matter of 
suppressing information in 0^^. Consider as a matrix with indices k, m. 
Make a Hermitian square matrix that sums over the information in index m: 

00t = ^ |e" > (A")2 < e"|; 

a 

< k\(j)(j)1\k' > = ^<7/.(/)|j'=|TO><m|j'='|V^(0>, 

m 

= K^^' . (13) 

The M-frames disappear in the sum due to unitarity. This is a rather magical 
version of unitarity, as the coset with — / < m < Z indices only just spans 
the three-dimensional subspace a = 1 ... 3 it needs to span. The diagonal frame 
of K coincides with the diagonal left-frame of </)^, namely the frame defined by 
the e's. Meanwhile the weights of K are the weights of (/) squared. 

We have shown, then, that a potential objection to an axis obtained from 
one particular method quadratic in the data is not a worry, because it is the 
same axis inherent linearly in the data. The objection remains true for other 
methods. 

2.3 Our Calculation 

We made our calculation (Table ^ for the quadrupole and octupole by making 
the spin-2 and spin-3 angular momentum matrices and diagonalizing by svd. 
The results are compatible with dOZTH with small differences. Wc find frame 
axes with largest eigenvalues: 

e(i)(Z 2) = (0.112, 0.506, -0.855); spin-2 

e(i)(Z = 3) = (0.246, 0.399, -0.883). spin - 3 (14) 

Our sign conventions are fixed by A" > 0, which is standard. The eigenvalues 
and complete frames are listed in Tabled We checked our result by comparing 
them to a direct diagonalization of the traceless symmetric Cartesian tensor for 
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spin 


A" 


ei 


62 


63 




25.95 


0.1115 


0.5049 


-0.8559 




17.70 


0.3874 


0.771 


0.5053 




8.24 


0.9151 


-0.3880 


-0.1095 


1=3 


47.35 


0.2462 


0.3992 


-0.8831 




21.94 


- 0.9383 


0.3262 


-0.1141 




21.24 


0.2425 


0.8568 


0.4550 



Table 1: Cosmic Axial Coincidence of quadrupole and octupole orientation. 
Top Group: Frame matrix elements from the quadrupole for e^(/ = 2) in rows 
a = 1..3 and columns fc=1...3. Bottom Group: Same for the octupole e'^{l = 3). 
The octupole has a nearly perfect degeneracy, allowing the lower two rows to 
be freely rotated about the e^{l = 3) axis. In that case the two frames nearly 
coincide. 

spin-2 , and by a third method using Clebsch-Gordan coefhcients (see below). 
We believe the small differences are not significant: n{l — 2)-e'^^^(Z = 2)=0.9996. 
There will be small differences going from the pixels to the moments, and it is 
impossible to obtain very high accuracy for an eigenvector by a variational 
method due to the object function becoming stationary. 

What about the eigenvalues? Consulting Tabled the largest eigenvalue is 
indeed large, and apparently inconsistent with an uncorrelated isotropic pro- 
posal. We will discuss quantifying such features shortly. The octupole's two 
"small" eigenvales A'^^^ A*^^^ are nearly degenerate. The frame of the octupole 
e^(/ = 3) is far from isotropic but can be aligned with the other frames with 
very small changes in the inputs. 

2.4 Other Methods 

We gave thought to the question of uniqueness, because it is tedious to analyze 
what other axial quantities might be made using higher powers of variables 
beyond the linear and quadratic constructions discussed so far. 

Under the rotation group, the product of two or more representations |^, m > 
m' > is reduced to a sum of tensors of rank \ J, M > by the Clebsch-Gordan 
series or "clebsches". So if we want a quadrupole in order to make axes, and 
it must come from from the product of two octupoles (say), then we need the 
3 (8) 3 ^ 2 clebsches and it can be worked out. We must distinguish the maps 
of \l, m > m' > and \l, m > (g) <', m'|, of course. 

2.4.1 Products of ^ (g) T 

The clebsches for reaZ symmetric projections l^l* 1 vanish, just as the expec- 
tation < J >= indicates. In fact J is just the operator of clebsches to make 
spin-1 from such products. We do not seek a vector, but a symmetric rank-2 
tensor to make axes. Evidently C'"'"' (2) =< j'^j'"' - \ j^5^^ > is the appropri- 
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ate map from 1(^1* ^ 2, as just discussed. The unit operator \ J'^8kk' added to 
make a traceless irreducible representation does not change the eigenvectors, but 
does remove one independent eigenvalue. The connection with clebsches makes 
the procedure of dOZTH less miraculous and shows that the svd frames, then, 
coincide with a Clebsh-Gordan series of squaring the data, projecting down into 
spin 2, and finding the preferred frame of the product. 

Although the result is trivial we checked our calculation of n2,'^3 a third 
way by generating it by the CG series for I ®l* . 

2.4.2 Products oil® I 

What about the clebsch projections of products without complex conjugation? 
One finds that l®l — > 1 and l^l 2, and so on, are non-zero, and perfectly valid 
complex quantities. If one's experience is limited to neutrino physics, these are 
like "Majorana" masses which take the non-conjugated products of the Lorentz 
group and make a scalar. From the svd decomposition one can gain insight on 
what they mean. Form the outer product 

C-^™' = E ete^A'^A'^ul-uJ (15) 

a, 13 

Now a clebsch series to get spin- J will project this by contracting with a matrix 
^mm' ■ '^^^ overall form resulting is 

QwC'/'I"'=e^W^fef,. (16) 

This is far from a diagonal form as there is no "dagger" on either symbol e. 
The function of W"^ is to mix the information in the wj„ not transforming like 
vectors until spin-J comes out of the mix. 

While spin-2 can be made, it is complex, and it is hard to get excited about 
an TO = 1 (say) term made from the product of / = 97, to = 93 and I — 97, to = 
94, TO = 92 that never existed before the data was squared. We do not pursue 
this or further products, given the invertibility of 0^ (Eq. that tells us it is 
complete. 

3 Cosmic Axial Coincidences 

We constructed symmetry methods to pursue a higher goal than classification. 
We turn to several striking coincidences among the CMB.2 and other axial 
quantities, all obtained from radiation propagating on cosmological distances. 
The coincidences pose a significant puzzle. 

First, the quadrupole and octupole axes coincide closely with each other, 
and with a well-known third axis e^ipoi^ pointed toward Virgo. The relations 
are: 

^dipole ' ^quad 0.980 
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6(jMod ' ^octo — 0-985 
^octo ' ^dipole 0.939. 



(17) 



Compare these to the default symmetry-based hypothesis that each vector 
is uncorrelated with the others and distributed isotropically. This happens to 
be the basic Big Bang (bBB) null, and no further assumptions are needed to 
employ it. The coincidence of three vectors in nearly the same direction is not 
at all expected. 

3.1 Coincidence with Independent Radio Frequency Ob- 
servables 

In Ref. [3] an axis of anisotropy for radiation propagating on cosmological scales 
was reported: in galactic coordinates the axis is found to be 

sjn = (-0.0927,0.4616,-0.8822). 

This axis coincides with the direction of Virgo. It is very remarkable that this 
axis coincides with axial parameters from the CMB. 

The result came from a test of propagation isotropy. The paper examined 
statistics of offset angles of radio galaxy symmetry axes relative to their aver- 
age polarization angles for 332 sources[3]. The polarization offsets are variables 
understood to be independent of intervening magnetic field. Faraday rotation 
measures RM are removed source-by-source in linear fits to wavelength depen- 
dence. Further discussion is given in the Section on Mechanisms. Now there 
are four coincident axes. 

3.2 Coincidence with Optical Frequency Observable 

In Refs. 0] Hutsemekers observed that optical polarization data from cosmo- 
logically distant and widely separated quasars indicate an improbable degree 
of coherence. The alignment of polarizations on a curved manifold (the sky) 
is tricky to assess and certain statistics invariant under coordinate changes are 
needed. In Ref. 5 invariant statistics were constructed using parallel-transport. 
The facts of unusual coherence were confirmed, and moreover a significant clus- 
tering of polarization coherence in large patches in the sky was observed. 
The axis of correlation obtained from a fit is 

SHut = (0.0344,0.4844,-0.8741) , 

in galactic coordinates. The axis given above has been obtained by imposing 
the cut on the degree of polarization p < 2 %. This cut gave the largest optical 
polarization alignment effect j5,. The axis obtained by the full data agrees 
within errors. This axis coincides with the direction of Virgo. 
Now there are five coincident axes. 
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3.3 CMB Polarization Observables 



COBE led the way in measuring moments of the temperature map beyond the 
dipole. WMAP |S] observations yield highly detailed information on CMB po- 
larizations for the first time. Among the features not predicted by the bBB is a 
highly significant cross-correlation of the temperature and degree of polarization 
(Stoke's Q) across very large angular scales. 

WMAP uses an invariant polarization statistic equivalent to parallel transport [S]. 
The observed correlations are inconsistent with the bBB predictions of uncorre- 
lated isotropy on such scales. A subsequent re-ionizing phase has been proposed 
to explain them. 

We remarked earlier that the conventional bBB has been ruled out. Insertion 
of a new hypothesis to explain the polarization correlation justifies that remark. 
Meanwhile we see no reason, even under the burden of the new hypothesis, for 
the polarization quantities to have any preferred axial character. 

For some reason WMAP has not released polarization data, as far as we 
know. We predict that WMAP may see an increased degree of polarization in 
the direction of Virgo, or diametrically opposite. One cannot deduce this from 
the correlation itself. We simply see a pattern here. Our prediction is empirical 
and based on the fact that the CMB already points to Virgo in the dipole, 
quadrupole, and octupole axes. Correlation with Virgo would simply be more 
of the Virgo Alignment. 

3.3.1 Directionality 

Given the opportunity to speculate, the pattern of data also indicates to us 
that CMB-derived polarizations should be aligned towards Virgo, along the 
plane found by Hutsemekers. Like the others, this deduction does not follow 
from the correlations published. It follows on the basis of inspecting patterns 
in data without imposing too many prior assumptions. 

3.4 Quantification with Bayesian Priors 

Here we consider the statistical significance of some coincidences. 
The null is uncorrelated isotropy. The usual test of the null is 



The probability of two axes Ca, &b coinciding from the null to a degree equal 
or better than the data is the area of the largest spherical cap touching them 
divided by 27r: 



P{null I data) 



P{data I null)P{null) 



(18) 



P{data) 



P{eA ■ es I null) = 1 - |eA • es]. 



(19) 



The probability of three coincidences is the product: 



P{eA ■ BB, CB ■ ec\ null) = ( 1 - |eA • esl )( 1 - |es • ec\) )■ 



(20) 
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P(null) 



Figure 2: Updated probabilities disfavoring the conventional big-bang isotropic 
null hypothesis. The solid curve shows the calculated probability of the hypoth- 
esis logio{P{null I Saxes) ) after seeing the coincidence of dipole, quadrupole, 
and octupole axes (only) as a function of prior prejudice P{null). The dashed 
curve shows logw{l — P{null)), an estimator of the probability of the alternate 
hypothesis. 
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dipolc quad octo JR 



Hutsemekers 



dipolc 
quad 
octo 
JR 



0.020 0.061 0.042 
0.015 0.023 
0.059 



0.024 
0.004 
0.026 
0.008 



Table 2: Probabilities as P-values P, for the axis coincidence listed to equal or 
exceed that seen in the data under an uncorrelated isotropic distribution. By 
ordinary convention 1 — P is the confidence level that the uncorrelated isotropic 
hypothesis is false. 

We summarize the numerical value of pairwise probabilities in Tabled 

Continuing, we need P{null), which comes from outside, and is difficult to 
settle. There may be an impulse to revert to large / correlations and declare 
P{null) = 1. But large I > 100 correlations made in a manifestly isotropic 
way add no independent information to the problem at hand. The things com- 
monly predicted, namely magnitudes of the low I moments, don't agree with 
the standard null. We find there is no good value for P{null). 

For argument let P{null) — 0.999 simply because we can say it is so. We 
need an alternate hypothesis P{axes) to proceed. This will be the hypothesis 
that axes are well-correlated. It may appear retrospective, but follows from 
the non-CMB data (Section 3). Quibbling is pointless because in a two-option 
contest P{axes) = 1 — P{null) = 0.001 << 1 is very unlikely anyway. By 
definition P(eA ■ gb, gb ■ ec) ~ 1 if axes tend to be correlated. So long as 
P(axes) is small enough the details of P{axes) will not change the evaluation 
of the data. Now calculate by Bayes Theorem: after seeing the coincidence of 
the dipole, quadrupole, and octupole axes, the updated probability of the null: 



The null survives because we used such a strong prejudice favoring the bBB. 
Yet if one had expressed less confidence, say P(null) = 0.990, the result would 
have changed, yielding P{null\3axes) 0.029, and ruling out the null. As 
a crisp summary, the null is ruled out with 97.1% confidence level, or indeed 
much more, unless one's prior prejudice favoring it is 99.0% or higher. For 
convenience we show parametrically the data-driven confidence level in the null 
and the alternative hypothesis, as a function of initial prejudice (Fig. O. 

3.4.1 Other Arguments 

It is interesting how a coincidence can be explained away by a Bayesian re- 
arrangement of logic. Assuming from the boost interpretation that the dipole 



P{null |3 axes) 



P(3 axes \ null)P{null) 



P(3 axes \ null)P{null) + P(3 axes \ cor)P{cor) ' 
0.02 X 0.015 X 0.999 



0.02 X 0.015 X 0.999 + 1 x 0.001 

= 0.23. 



(21) 
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axis "can be of no cosmological significance" (goes the argument), we remove 
one vector freedom in the product of three possible probabihties. Then only two 
things coincide, and the astonishing improbability of the coincidence is reduced 
to a merely alarming value. 

Is this a swindle? It depends. Nobody can argue with a Bayesian argument 
because the proponent who defines the priors may define their own priors. The 
simpler and more direct expression of prior prejudice will hold that all the low 
/ moments have no (or little) cosmological significance. In that case, nobody 
should not have even looked at (let alone paid millions for) the data! 

There is a Bayesian response: let P{junk) be the probability that \ow-l is 
junk. "Cosmic variance" , in the assertion that the fluctuations of a few degrees 
of freedom tend to equal their value, is for us on the logical track proportional 
to P{junk). In that case the comparison of three axes would not be meaningful, 
except to update prior beliefs in P{junk). 

We are not ready to assume that the bulk of the data in a large-scale sense 
is just extremely expensive junk. Given the cost and the potential importance 
of the decision we feel it would be highly irresponsible to do so. The logical 
chain proportional to 1 — P{junk) is the chain capable of testing the null. The 
smallncss of 1 — P{junk) is the entrance price for asking challenging questions. 
One is not supposed to assess the entrance fee a second time when something 
surprising is actually found. 

Under the original hypothesis, and not adding further qualifications, the 
conditional probability of coincidence P{ni, n2,fi3,\ bBB) of three axes is an 
extremely unlikely occurence. 

4 Mechanisms 

The optical properties of cosmological propagation are not directly connected 
to the Big Bang by general principles and need not have significant other con- 
sequences. 

It is also not necessary to have a mechanism to recognize regularities in data. 
Proposing mechanisms is helpful toward suggesting interpretations and further 
tests. 

Here we review information on the non-CMB axes and then turn to potential 
mechanisms. 

4.1 Faraday Rotation Data 

There is a long history of observation of anisotropy at GHz frequencies. Contro- 
versies have been generated by discordant signals in statistics initially supposed 
to represent the same phenomena. Over a period of years, the discrepancies 
were traced to different symmetries in the quantities being reported. There is a 
consistent signal in odd-parity statistics sensitive to a direction of polarization 
rotation. 
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X = sin(2 beta) 



• 


X > 0.66 


• 


0.66>x>0.33 


• 


0.33 >x>0 




> X > -0.33 


• 


-0.33 > X > -0.66 


• 


X < -0.66 




Figure 3: Covariant polarization offset statistics x = sin{2f3), where /? is the 
offset angle between galaxy symmetry axes and polarization planes, for 265 radio 
galaxy sources mentioned in the text. The prefered axis {I — 101.36° ± 15°, b = 
—61.91° ± 15°) is shown as a black cross. Axial clustering is readily visible to 
the eye. 
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Birch^ observed patterns in data empirically. Kendall and Young^j con- 
firmed Birch with a likelihood analysis. Bietenholz and Kronberg^^j verified 
Birch's results, and confirmed a statistically significant effect in a larger data 
set. They dismissed the claim by turning to a statistic of even parity symme- 
try developed by the statisticians Jupp and Mardia^]. Nodland and Ralston 
[T^ reported a redshift-dependent odrf-parity correlation in a set of 160 sources 
with known redshifts. These data indicated correlations of offsets with an axis 
oriented^ toward Virgo. Carrol and Field J3] dismissed the signal using RMS 
angle measures that are coordinate-dependent^. Eisenstein and Bunn re- 
plotted the data on the basis of even parity and dismissed the signal by eye 
| 15II16 |. Loredo et al|17| found no signal in an invariant statistic of even parity. 

A symmetry analysis^J showed that the Jupp-Mardia {JM) statistics come 
in even and odd-parity classes. Ref. reported axial correlations obtained from 
a larger set of 332 radio galaxies. Integrating over redshifts made statistics more 
robust and increased the size of available data. Both coordinate-free odd-parity 
JM statistics and likelihood analysis showed a strong signal. The case of highest 
significance, about 3.5 cr, was found using an RM cut parameter.^ The cut is 
{RAI — RM) > 6 rad/rn^, where RM is the mean value of RM in the sample. 
A related study JH| yield axes consistent with the axis of Eq. 13. II within errors. 
In fig. 01 we show the scatter plot of sin(2/3) for the 265 source data sample 
considered in Ref. after imposing the cut {RM — RM) > 6 rad/m? . Here 
(3 is the angle between the observation polarization angle, after taking out the 
effect of the Faraday Rotation, and the galaxy orientation angle. The figure 
shows the averaged value of the variable x = sin(2/3) at each position, where the 
average is taken over the nearest neighbours which satisfy fi-f2> 0.95. Here 
fi and 7^2 refer to the angular positions of the two sources. 

The mechanism for these correlations remains unknown. The population in 
RM is peculiar, with a tall central spike and long tails. The population and 
correlations might suggest an artifact of the data-fitting procedure, which used 
the faulty concept of RMS angle. However re-fitting the original data to obtain 
RM and offset values by an invariant statistical method ^HI did not change the 
conclusions significantly. 

4.1.1 Plasma Explanations? 

Perhaps Faraday rotation results could be explained away by the ineluctable 
peculiarities of plasma electrodynamics. One study of different statistics taken 
at higher frequencies saw no signal in a selection of polarization and jet axes. 

1201 

Faraday rotation depends on the longitudinal component _B|| of the mag- 
netic field, which must certainly fluctuate over cosmological distances. However 
adiabatic propagation models traditionally used break down at field direction 
reversals |2] B\\ 0. This conventional physics was not included in the astro- 
nomical RM analysis, and flips polarizations. It might be a candidate for ex- 

•^The direction diametrically opposite to Virgo was reported, but spfn has no sign. 
*The statistical effect of choosing an RM cut parameter is of course included. 
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planation. Large angular scale correlations in galactic or cosmological magnetic 
fields, however, would be required for this mechanism to explain the Faraday 
rotation data. Galactic fields were always a concern for RM measurements and 
were thought to be understood. We continue along these lines in the subsection 
on Foreground Mechanisms below. 

Cosmological magnetic fields of dipole, quadrupole or octupole character 
have always formed a highly virulent threat to Big Bang dynamics. Little re- 
liable information exists about cosmological magnetic fields j22j, and data at 
frequencies dominated by Faraday rotation poses inherent difficulties of inter- 
pretation. 

The CMB data, with frequency much larger than Faraday rotation measure- 
ments, are thought to be comparatively free of such effects. Compensations are 
made for the effects of the ionosphere and galaxy. Under the standard hypothe- 
ses, the CMB measurements form an independent observable which should not 
be related to Faraday rotation data. 

4.2 Optical Data 

The interpretation of Hutsemekers result and subsequent verification is compli- 
cated by the possible processes affecting optical polarization along the trajec- 
tory. The primary model to change the degree of polarization is preferential 
extinction by dust. However any correlation of dust across vast cosmological 
distances would be as alarming as a correlation of magnetic fields or the CMB 
radio spectrum. 

Peculiarly, the data fraction of smallest polarizations has the highest sig- 
nificance. In Fig. 2]we show the distribution of correlation statistic x — S^, 
defined in Ref. , for the 146 source data sample after imposing the cut p < 2 
%. The statistic is the mean value of the dispersion over the entire sam- 
ple. A measure of the dispersion at each point is obtained by maximizing the 
function, 

d^{0) = — y cos [29 ~2{ek + Afe^,)] , 

k=l 

with respect to the angle 9. Here {6k + ^k^i) is the polarization at the position 
k after parallel transport to position i along the great circle joining these two 
points. The symbol is the number of nearest neighbours. Fig. Pluses Uy — 28, 
a value giving the strongest signal of alignment. 

Although thought to be under control, we cannot rule out foreground effects 
here and refer to the subsection on Foreground Mechanisms below. 

4.3 Mechanisms 

We considered the following mechanisms: 

• Foreground Mechanisms: All of the data cited has potential weaknesses 
in depending on models of the foreground, and in particular, the effects of our 
Galaxy. 
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Figure 4: Optical polarization alignment observed by Hutsemekers 0] in (a) 
equatorial and (b) galactic coordinates. Here x refers to the correlation statistic, 
defined in text, with the number of nearest neighbours chosen to be 28. The data 
within the galactic plane was deleted 4 in order to minimize the effect of galactic 
extinction. Effects of population are taken into account in the construction of 
the prefered axis, shown as a black cross at I = 85.94°, b = —60.94°. Axial 
clustering is readily visible to the eye. 
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A great deal of ingenious technical effort has gone into taking this into ac- 
count without inducing biases. It is unlikely that statistical procedural errors 
exist. Nevertheless there may be systematic physical processes not taken into ac- 
count in current CMB foreground subtraction, which cause an alignment with 
features of our Galaxy. In the spirit of current CMB fits, there is reason to 
believe that free adjustment of many galactic magnetic field and optical extinc- 
tion parameters might explain not only the CMB data but also the non-CMB 
correlations. 

It is somewhat hard to see how substantial foreground corrections that might 
make axial alignment correlations go away would not change CMB data and 
preserve its current interpretation as a pristine record of the early universe. 

• Metric Mechanisms: The propagation of polarizations is an old subject in 
General Relativity In the freely-falling frame where the metric is triv- 
ial, general covariance says that the two polarizations propagate equally in 
Maxwell's theory, and will preserve the state and direction of polarization. This 
is due to a symmetry of duality that the metric theory happens to preserve. 
Large polarization effects would need large curvature effects that are not ob- 
served. Given present information and beliefs about the large-scale flatness of 
the Universe, wc do not know of a way for the large effects in the polarization 
data cited to be explained. 

• Magnetic Fields in the Early Universe: The scaling laws of magnetic fields 
are such that their imprint on the present Universe can be hard to trace. How- 
ever large-scale magnetic fields clearly could affect the CMB. We cannot review 
all the work done in this regard: It is simply hard to see how magnetic fields 
at small redshift could conspire in conventional physics to explain the Faraday 
offset anisotropics, and it seems impossible to explain the Hutsemekers data. 

• Dark Energy: If dark energy is due to a scalar^ field there is a simple generic 
mechanism to produce surprising optical effects that do not necessarily disturb 
the other predictions of bBB models. Since we have not seen this discussed 
elsewhere, we summarize it in broad terms, referring the reader to a companion 
paper 1231 . 

During the process of expansion and decoupling, the photon field changes 
dramatically over a wider range than any other particle. The change occurs in 
the photon mass. The photon is massless en vacuo due to the local nature of 
the theory and gauge invariance. The photon is not massless in general, but 
propagates with a mass given by the "plasma frequency" — m^c'^. The 
usual non-relativistic formula for the plasma frequency is 

'^P 47rme ' 

where — e is the electron's charge, n is the density of light (electron) carriers, 
and rrig is the electron mass. The dispersion relation en plasmo is 

= my. (22) 

^As common in cosmology we include pseudoscalars with sealars. 



20 



Any scalar field (p coupling to the electromagnetic field strength F^"^ will 
mix resonantly in a background magnetic field when the scalar mass is same as 
the photon mass. A typical Action is 

S = j d^x^ j^-l^'^-F^^+g^e^^^^F'^-^"/'^ , (23) 

which suffices to break duality symmetry. Parameter g is a dimensionful cou- 
pling constant that can be of exceedingly small size. The physical consequences 
of the mixing of this scalar field with photons has been studied by several au- 
thors Since the plasma frequency and the scalar mass are independent, and 
of different scaling characteristics in general, it is hard to avoid scale-crossing. 
In particular, driving the photon mass to zero crosses every conceivable scale at 
least once, to say nothing about re-ionization phases. Once duality symmetry 
is broken, polarization can be spontaneously generated and also rotated under 
propagation. 

No fundamental coupling is required because radiative corrections will gen- 
erate all couplings allowed by symmetry. The existence of a scalar field -that is, 
dark energy - is itself sufficient to guarantee that the optical properties of prop- 
agation over vast distances and times will not be those of the usual theory. Of 
course, the organization of dark energy over huge scales is something that can- 
not be determined in advance. Models such as quintessence hardly encompass 
the possibilities. 

The observation of propagation anomalies is in fact a highly conservative pre- 
diction of the current framework of cosmology, without requiring new hypothe- 
ses. While constantly cited, the data has ruled out Einstein's early cosmological 
constant, the dark action 

Sdark = J d'^X^A. 

Instead it is necessary to have A — > A{x, t). Now it would violate basic con- 
sistency requirements, such as causality, if what is currently meant by A{x, t) 
is not associated with the modes of a field. We believe that what is meant by 
A{x, t) is at least as dynamical as a field even more than we believe in A(a;, t) 
Unfortunately the current state of development of theory cannot relate any 
particular field to any measurement to A(z, t). One wonders why it is so im- 
portant. Meanwhile ordinary perturbative methods can relate any field through 
radiative corrections to an action such as Eq. |2S1 Thus propagation anomalies 
may be an effective probe of dark energy, as compared to the rather indirect 
coupling to gravity via an energy-momentum tensor. We reiterate that light 
propagation may yield a productive way to look for the effects of dark energy. 

• Dark Energy Plus Magnetic Fields 

Given that a resonant interaction inevitably occurs between dark energy 
and electrodynamics, we have to admit ignorance on the potential outcomes in 
terms of large scale magnetic fields. This is an extremely interesting dynamical 
problem. The interactions of magnetic fields are not expected to retain isotropy 
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in all features of an expanding Universe. Father study along these directions 
seems justified. 

5 Concluding Remarks 

There is nothing to fault in the beautiful technology bringing moments of the 
CMB well above I > 100. Yet the small moments are problematic. The Virgo 
alignment appears to contradict the standard picture. 

We believe that the tests possible in the near future and longer term hinge 
on polarization quantities. Perhaps the best way to search for dark energy is 
with light. WMAP may be poised to make a major discovery that goes far 
beyond the verification of the models of previous generations. 

6 Appendix: Conventions 

Here we collect details of conventions and procedures. 
For spin I we define operators 

Jz — {—I, —I + 1, ■■■I — 1, I) (diagonal), 

and we enter multipole moments as column vectors in the same convention. We 
define 

J, = U{0, -7T/2, 0)J,[/t(o, -7r/2, 0) , 

where the matrix elements of U{0, /3, 0)mm' = dl^m'W Wigner d- 

matrices with the standard phase conventions. Finally, 

Jy = i [Jz^ Jx ]• 

The Hutsemekers axis is obtained as follows. Let fij be the unit vector 

coordinates of the ith source on the sky in galactic coordinates. Construct a 
two point comparison between the i and j polarization, given by 

C{ni,nj) = cos[29i - 2{0j + Aj^i)] , 

where {9j + Aj^i) is the polarization angle obtained after parallel transport 
of the polarization plane at position j to the position i along the great circle 
connecting these two points. Construct the correlation tensor T'^'^ , 

By construction T — > in a population of uncorrelated polarizations and posi- 
tions. Make the matrix 

s = m-i/2tm-^/2 
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where M is the matrix of sky locations 



By construction S ^< p > 1 in an isotropic sample with a net polarization 
parameter < p >, so that S* — > in the null. The Hutsemekers axis is the 
eigenvector of S with the largest eigenvalue. 
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